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We compute the isospin-asymmetry dependence of microscopic optical model potentials from 
realistic chiral two- and three-body interactions over a range of resolution scales A ~ 400 — 500 MeV. 

We show that at moderate projectile energies, ^inv = 110 — 200 MeV, the real isovector part of the 
optical potential changes sign, a phenomenon referred to as isospin inversion. We also extract the 
strength and energy dependence of the imaginary isovector optical potential and find no evidence 
for an analogous phenomenon over the range of energies, E < 200 MeV, considered in the present 
work. Finally, we compute for the first time the leading corrections to the Lane parametrization for 
the isospin-asymmetry dependence of the optical potential and observe an enhanced importance at 
low scattering energies. 


Introduction - The structure and dynamics of neutron- 
rich nuclei are key inputs for modeling neutron stars, 
core-collapse supernovae and r-process nucleosynthesis 
[1-8]. Elucidating the properties of highly isospin- 
asymmetric nuclear matter is therefore a priority in low- 
energy nuclear science research and a major motivation 
for the development of next-generation radioactive ion 
beam (RIB) facilities. Microscopic many-body meth¬ 
ods [9-11] with chiral two- and three-body forces have 
been successful in describing the bound-state properties 
of neutron-rich matter. Complementary and consistent 
nuclear reaction models are under development [12-14], 
and of these, global optical potentials aim to address the 
broadest theory needs for interpreting RIB scattering ex¬ 
periments and simulating r-process nucleosynthesis. In 
fact, current modeling of the strong r-process favors a 
cold scenario in binary neutron star mergers [6-8, 15- 
17], where mass transfer to highly neutron-rich isotopes 
occurs and freeze-out is achieved more rapidly, which en¬ 
hances the importance of radiative neutron capture pro¬ 
cesses in determining the final abundance pattern of r- 
process elements [18]. 

Although global phenomenological optical potentials 
[19-21] are well constrained close to the valley of nu¬ 
clear stability, their predictive power for reactions in¬ 
volving exotic neutron-rich isotopes is not well under¬ 
stood. Elastic scattering data has been used in the past 
to parametrize local optical potentials in specific regions 
of the nuclear chart off stability, however, the most exotic 
and low-intensity radioactive ion beams require thick- 
target experiments that provide quality inelastic data 
only [22]. This motivates the need for accurate micro¬ 
scopic optical model potentials and investigations of their 
energy and isospin-asymmetry, 6np = {N — Z^jA, de¬ 
pendence. Identifying energy regimes in which the lead¬ 
ing linear term is dominant can then be valuable for 
extrapolating existing phenomenological potentials away 
from stability. 

The aim of the present work is to employ high-precision 
chiral two- and three-nucleon forces to study the real and 


imaginary volume components of the nucleon-nucleus op¬ 
tical model potential far from the valley of stability. The 
dependence on the isospin asymmetry of the target nu¬ 
cleus is traditionally taken to be linear and isovector in 
character, a parametrization known as the Lane form 
[23]. The isoscalar components of the optical potential 
are then independent of d^p^ We revisit these assump¬ 
tions and find that subleading terms in d^p (which are 
isoscalar and isovector for even and odd powers of dnp^ 
respectively) can be important for highly neutron-rich 
nuclei and particularly at the low energies most relevant 
for nuclear astrophysical phenomena. We study the en¬ 
ergy dependence of these terms up to E 200 MeV in 
anticipation of future experimental investigations of ex¬ 
otic isotope reactions at RIB facilities. 

A phenomenon of particular interest is isospin inver¬ 
sion, whereby the real isovector optical potential is ex¬ 
pected to change sign from positive at low energies to 
negative at higher energies. In the vicinity of isospin 
inversion, subleading terms proportional to may be¬ 
come relevant and probe novel isospin physics. The in¬ 
terplay between intermediate-range attractive contribu¬ 
tions to the nucleon-nucleon interaction and short-range 
repulsive contributions can give rise to a change in sign 
of the isoscalar optical potential at projectile energies 
greater than E ^ 250 MeV, which has been observed 
in calculations with microscopic nucleon-nucleon poten¬ 
tials [24]. The isovector contribution to the optical po¬ 
tential, arising from 7r-meson exchange and p-meson ex¬ 
change in traditional one-boson exchange models, has a 
stronger relative energy dependence. In the past, semi- 
microscopic and microscopic optical potentials have been 
constructed from mean field theory [25-27] and realis¬ 
tic nucleon-nucleon interactions [28-33] respectively, and 
there is significant disagreement regarding the energy de¬ 
pendence of the isovector components. To date there are 
therefore no strong constraints on the kinematic transi¬ 
tion region associated with isospin inversion. 

In microscopic many-body theory the nucleon-nucleus 
optical model potential is identified with the on-shell nu- 
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cleon self-energy S(r, r';^) [34]. In the present study 
we compute the nucleon self-energy at second order in 
many-body perturbation theory employing as a start¬ 
ing point high-precision nuclear interactions derived from 
chiral effective field theory [35, 36]. Chiral nuclear poten¬ 
tials with momentum-space cutoffs A < 450 MeV [37-39] 
exhibit very good perturbative behavior (comparable to 
renormalization-group evolved potentials [40-42]) and we 
also consider a potential employing a 500 MeV cutoff [36] 
that is used to give a conservative estimate of theoretical 
uncertainties associated with nonperturbative dynamics 
and variations in the resolution scale. Previous work has 
shown that chiral low-momentum potentials provide a 
good description of the symmetric nuclear matter sat¬ 
uration energy and density [39], the incompressibility 
and isospin-asymmetry energy [43], and the critical end¬ 
point of the liquid-gas phase transition [44, 45]. The 
present approach to nuclear scattering is therefore con¬ 
sistent with nontrivial constraints from nuclear structure. 

Iso spin-asymmetry dependence of optical potentials - 
In the optical model for nucleon-nucleus scattering, the 
complicated many-body problem associated with mul¬ 
tiple scattering through two- and three-body forces is 
replaced by an (elastic-scattering) equivalent complex¬ 
valued single-particle potential: 

V(f, f E) = U{r, f E) + iW{r, r';E), (1) 

which in general is both non-local and energy-dependent. 
The imaginary part in Eq. (1) accounts for the presence 
of open inelastic scattering channels. Phenomenological 
optical potentials are often taken to be local and energy- 
dependent. The real and imaginary parts contain vol¬ 
ume components proportional to Woods-Saxon densities 

fj{r) = + 

U{r-,E) = -U{E)fr{r), W{r-, E) =-W{E)fi{r), ( 2 ) 

where the parameters Uo{E),Wo{E),Rr,i and ar,i vary 
smoothly with the mass number A of the target nucleus 
and the projectile energy E. 

Recently chiral two- and three-nucleon forces have 
been used to compute the real and imaginary volume 
components of the optical potential for isospin-symmetric 
nuclear systems [12, 13]. Although the strength and 
energy dependence of the real component was found 
to be in good agreement with modern phenomenologi¬ 
cal parametrizations [21], the absorptive strength of the 
imaginary part from microscopic nuclear potentials was 
about a factor of two larger than that from phenomenol¬ 
ogy. We note that at low energies the phenomenological 
imaginary part is surface peaked and vanishes in the infi¬ 
nite matter limit. Microscopic many-body theory in the 
local density approximation attributes the surface imag¬ 
inary part to the nonlinear density dependence of the 
imaginary volume component. In Fig. 1 we show the real 
and imaginary parts of the optical potential at nuclear 
matter saturation density po = 0.16 fm“^ from chiral nu¬ 
clear potentials compared to the global fit in Ref. [21]. 



FIG. 1: (color online) Energy dependence of the real and 
imaginary parts of the microscopic optical model potential 
from chiral two- and three-body forces for symmetric nuclear 
matter at saturation density po- Shown for comparison are 
the global phenomenological potentials of Ref. [19-21]. 

The phenomenological “Koning” bands are obtained by 
varying the mass number over the range A = 50 — 150 and 
should not be interpreted as an uncertainty. On the other 
hand, the error band associated with the microscopic cal¬ 
culation comes from varying the resolution scale over the 
range A = 414 — 500 MeV. In Fig. 1 we have included 
for comparison also global optical potential parametriza¬ 
tions [19, 20] valid at lower energies with associated un¬ 
certainty estimates. In contrast to the results reported 
in Ref. [12], the single-particle energies entering in the 
second-order perturbative calculation are computed self- 
consistently via e{q) = g^/2M + Re E(g, e(g)) rather than 
from the effective mass plus energy shift parametrization 
e{q) g^/2M* + A, which smears out the enhancement 
of the momentum-dependent effective mass at the Fermi 
surface [46]. From Fig. 1 we observe that the microscopic 
real volume component remains nearly linearly depen¬ 
dent on the incident energy beyond E = 100 MeV, in 
contrast to the Koning analysis in Ref. [21]. Neverthe¬ 
less, the two approaches are consistent within uncertain¬ 
ties over a wide range of energies. 

For scattering on neutron-rich nuclei, the dependence 
of the optical potential on the isospin asymmetry 6np is 
crucial. The standard Lane form 

U = Uo + ^Ui, (3) 

where r and T are the isospin operators for the projectile 
and target nucleus respectively, is widely used in both 
phenomenological and microscopic calculations. The 
Lane parametrization relates the elastic proton-nucleus, 
neutron-nucleus, and quasi-elastic charge-exchange pro¬ 
cesses. For elastic scattering the Lane form reduces to 
U = Uo — UjdnpTs, where rs = ±1 (for protons and neu¬ 
trons, respectively) is the isospin quantum number of the 
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FIG. 2: Diagrams contributing to the proton and neutron 
self energies Ep,n(^, Ct;;/c/, (5np) at first and second order in 
perturbation theory. The wavy line represents the antisym¬ 
metrized two-nucleon interaction V 2 N and the thick wavy lines 
on the second-order diagrams represent the sum of the free- 
space two-body force and the density-dependent NN interac¬ 
tion from Refs. [47, 48]. 
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incident nucleon. 


Here we consider a more general expansion of the 
isospin asymmetry dependence: 

U = Uo-Ui nSnp + UiiSlp + 0{6lp). (4) 

The Hartree-Fock contribution {q\kf,5np) from 

two-body forces, shown diagrammatically in Fig. 2(c), 
is obtained by generalizing the results of Ref. [12] and 
has the form 

= '^{qhissitti\V2N\qhissitti)ni, 

( 5 ) 

where V 2 N is the anti-symmetrized potential matrix ele¬ 
ment, rip^n = T — |^i|) is the occupation 

probability, and the sum is taken over the momentum 
/ii, spin si, and isospin ti of the intermediate hole state. 
In the present work we compute the exact dependence 
of {q; kf ^ 6np) on and extract the linear and 

quadratic terms numerically. 


The Hartree-Fock contributions from three-body forces 
are obtained by summing two particles over the filled 
states in the fermi sea: 


^(i),3Ar 


(^, kf , 


(6) 


1 

2 


^ 1 ^ 2 ; SS 1 S 2 ; ttit2\V3N\qhih2; ssiS2;ttit2)nin2. 
12 


The next-to-next-to-leading order (N^LO) chiral three- 
body force consists of three terms, whose diagrammatic 
contributions to the nucleon self-energy are shown in Fig. 
3. The chiral three-nucleon contact force is proportional 
to the low-energy constant ce- 


V, 


(ct) 


3N 


= E 


Ce 


2/^A, 


‘3 ’ 


( 7 ) 


FIG. 3: Gontributions to the Hartree-Fock single-particle po¬ 
tential from the N^LO chiral three-nucleon force. The large 
dots represent vertices proportional to the low-energy con¬ 
stants Cl, C3, C4, cp), cp;, and the short double-line indicates a 
medium insertion: —27rS{ko)0{k^'^ — |^|). 


where = 700 MeV and = 92.4 MeV. The one-pion 
exchange three-body force proportional to the low-energy 
constant cd has the form 


w(l7r) _ _ 

i^j^k 


QaCd 




8/^A^ + ml 


' Qj 'Ti ' 'Tj , 


(8) 


where qa = 1-29 and = 138 MeV. Finally, the two- 
pion exchange component with vertices proportional to 
Cl, 3,4 is given by 




( 277 ) 


3N 


= E 

i^j^k 


9a 


Ci • Qi CTj • Qj 


8f^iqi +ml){qj + ml) 


potl3 a 3 
ijk ’ 


( 9 ) 


where the isospin tensor is 

(-4cim^ + 203 * • 9 i)+C 4 e“^'^r^CTfe-(* x *). 

( 10 ) 

The low-energy constants 01 ^ 3^4 have been fitted (within 
the empirical uncertainties imposed by ttN scattering) 
to nucleon-nucleon scattering phase shifts [36], while the 
Ce and ce constants have been fitted to reproduce the 
binding energies of ^H and ^He as well as the /3-decay 
lifetime of ^H [39]. 

At the Hartree-Fock level, the first-order terms in 6np 
arising from contact three-nucleon forces proportional to 
the low-energy constants ce and ce are given by 


Ul{q,kf) = 


csk^f 


9ACDm%v3 


{2u - 2A (11) 


3(27r/^)4Ax 
— arctan 2u — arctan(i4 + x) — arctan(rt — x) 

+ — In 1 + 4u^) +---In 

Au 12x 1 {u — xY 


where kj^ = 37r^(pn + Pp)l‘^i ^ and u = kf 

The 27r-exchange Hartree diagrams proportional to ci^s 
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give rise to the isovector optical potential strength func¬ 
tion 


Ui{q,kf) 




|c3W + 


Cl — C 3 1 + (u + a;)^ 
2x 1 + (m — x)^ 


(C3 - 2ci)m 


[1 + (u + a;)2][l + {u- a;)2] 


( 12 ) 


while the 27r-exchange Fock diagrams proportional to 
Cl, 3,4 yield 


Ui{q,kf) = 


g\mlu 


6ci 


u) duH{x^ u) 


9(47r/^)4x2 

+H{u,u) dxH{u,x) + ( 2 c 4 — C 3 ) G{x,u) duG{x,u) 
+G{u,u) dxG{u,x) — 2(c3+C4) I{x,u) dul{x,u) 

+I{u,u)dxI{u,x) + [ 

^ Jo 

+(3c3 + 2c4)duG{^, u) dxG{^, x) + 2{3cz - C4) 
xdul{^,u)dxl{^,x) 

with auxiliary functions 


18ciduH{^,u)dxH{^,x) 

(13) 


AiUX o 

G{x,u) = —— — 3) + 4x[arctan(?i + x) 

o 

/ /2 2 ^ {u x\^ 

+ arctanf 14 — x) -h (X —u — 1) In-^ 

1 -h (i4 — xy 


(14) 


H{x^ u) 


X 


14(1 -h 14^ + x‘^) — 
[l -h (i4 — In 


^[l + (« + x)2] 

1 + (i4 -h xy 
1 -h (i4 — xy ’ 


(15) 


I{x,u) = —(8i4^+3x^) 
6 

—+ (1 — 314^)(1 -h 14^ 



(l + w2)2 + ^ 

1 + (l4 -h x) 
^ 

1 + (l4 — X) 


(l+^")3 

x‘^ 

I ’ (16) 


The terms second-order in Snp are relatively small, and 
the explicit expressions are given in the Appendix. 

Finally, we consider the second-order perturbative con¬ 
tribution from two and three-body forces, b^ 2 ^+ 3 Ar’ which 
is approximated by employing a density-dependent NN 
interaction constructed from Vsn as described in Refs. 
[47-49]. The second-order contributions, shown in Figs. 
2(a) and 2(b), are given by 


^(2q,),2A7/ j £ \ 

^4 [Q") ^ ^np) 

_ 1 \{piP3SiS3titz\ViN\qh2SS2tt2)\‘^ 

2 ^ w + £2 - ei - e 3 + i?? 

x{2Tr)^5{pi+pz - q- h2), 


711712743 


(17) 
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FIG. 4: (color online) Contributions to the real part of the 
proton and neutron optical model potentials from the n31o450 
chiral two- and three-body forces as a function of the momen¬ 
tum q and isospin asymmetry 3np at nuclear matter saturation 
density po- 


{q,UJ-kf,8np) 

_ 1 |(/llh35l53^1^3|C2'^^|gP2gS2^^2)P 

^ m W + £2 - £1 - £3 - 777 

X ( 27 r)^( 5 (hi +hz-q-P2)- 


771772773 


(18) 


where h* = 1 — 77^ and is the sum of the free- 
space nucleon-nucleon potential V 2 N and the density- 
dependent two-body force ^ 2 ^*^ [47, 48]. All single¬ 
particle energies are computed self-consistently: 

ep,n(g) = q^/2M + (19) 

+ I]W-3^(g) + ReS(2)^(9,e(g)). 


The expressions are computed for arbitrary isospin asym¬ 
metry, and the linear and quadratic terms in Snp are ex¬ 
tracted numerically. 

We show in Fig. 4 the various contributions to the 
real part of the single-particle potential for protons and 
neutrons in asymmetric matter at a density of po com¬ 
puted from the n31o450 chiral NN potential. We notice 
that when employing the n31o450 low-momentum chi¬ 
ral nuclear potential (and also the n31o414 potential), 
the second-order perturbative contribution is always less 
than the Hartree-Fock contribution, in contrast to the 
behavior observed in Ref. [12] using the n31o500 poten¬ 
tial. We find that three-nucleon forces in general enhance 
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FIG. 6: Energy dependence of the isovector imaginary optical 
potential at saturation density from chiral two- and three- 
body forces. Also shown are the subleading contributions 
to both the real and imaginary potentials. 


FIG. 5: (color online) Energy dependence of the isovector real 
optical model potential at saturation density from chiral ef¬ 
fective held theory. Shown for comparison are the predictions 
of other microscopic, semi-microscopic, and phenomenological 
models (see Ref. [33]). 


isospin inversion due to the repulsive character of three- 
body forces in homogeneous matter. 

The magnitude and energy dependence of the real 
isovector part of the optical potential are poorly con¬ 
strained by experiment. From Refs. [19-21, 25, 50-52] 
one hnds that the magnitude is expected to decrease with 
energy according to Uj = (28 ± 6) MeV — (0.15 ± 0.05)E, 
which we show as the empirical band in Fig. 5. The 
chiral effective held theory prediction shown in Fig. 5 
is consistent with the empirical constraints but has sig- 
nihcantly smaller uncertainties, typically of about 5-10 
MeV over a wide range of scattering energies. The re¬ 
gion for isospin inversion is predicted to lie in the range 
^inv = 155 ± 45 MeV, and the two low-momentum in¬ 
teractions alone would give a much narrower region of 
^inv = 120 ±10 MeV. The results from other microscopic 
many-body calculations are shown, including Brueckner- 
Hartree-Fock “BHF” [53] and Dirac-Brueckner-Hartree- 
Fock “DBHF” [33], as well as semi-microscopic mean held 
models: “RMF” [54], “Gogny” [55], “SKM*” [56], and 
“SkLya” [57]. We refer the reader to Ref. [33] for addi¬ 
tional details and analysis. 

In addition to variations in the resolution scale, also 
the order-by-order convergence [58, 59] in the chiral ex¬ 
pansion provides an estimate of the theoretical uncer¬ 
tainty and associated cutoff artifacts. For this purpose 
we have computed as well the optical potential from the 
NLO (next-to-leading order) and N2LO chiral potentials 
with cutoffs A = 450 and 500 MeV [58]. At thresh¬ 
old the uncertainties for both sets of cutoffs are similar: 
Uf^{0) = 32±6 MeV and Uf^{0) = 31±6 MeV. At the 
isospin inversion energy the uncertainties are slightly re¬ 
duced: Uf^{Einv) = 0 ± 5 MeV and Uf^{Einv) = 0 ± 4 
MeV. Accounting for these uncertainties would not qual¬ 


itatively alter the error bands shown in Fig. 5, except 
at low scattering energies. We also note that beyond 
an energy of E 200 MeV, signihcant artifacts were ob¬ 
served in the calculation of the optical potential from the 
n31o414 potential, suggesting a breakdown in the chiral 
effective field theory expansion. 

In Fig. 6 we show the subleading contribution Un 
to the real optical potential as a function of projectile 
energy. For low-energy scattering on neutron-rich tar¬ 
gets we can expect an isoscalar shift of roughly (15 — 
20MeV)(5^p, while for energies greater than E > 100 
MeV the quadratic term is consistent with zero. The lat¬ 
ter observation has the important consequence that ex¬ 
trapolations of phenomenological optical potentials into 
the neutron-rich region of the nuclear chart can be valid 
for energies beyond > 100 MeV. This is supported by 
a similar feature in the volume imaginary component, 
where in Fig. 6 we see that the quadratic S^p term Wu 
is consistent with zero for all scattering energies consid¬ 
ered. The Lane parametrization of the volume imaginary 
optical potential strength therefore provides an excellent 
approximation to the true isospin asymmetry dependence 
over a large range of energies. We note that since no 
isovector component for the volume imaginary contribu¬ 
tion could be extracted from the most recent analyses in 
Refs. [20, 21] due to the uncertainties in the scattering 
data at large energies, our result is a prediction that may 
be verified at RIB facilities. 

Finally, we consider the quality of fitting the isospin 
asymmetry dependence of the proton and neutron real 
optical potentials up to quadratic 6‘^p terms. In Fig. 7 
we show as a representative example the isospin asymme¬ 
try dependence of the optical potential for a scattering 
energy of 50 MeV employing the n31o450 chiral nuclear 
potential. For values of the isospin asymmetry up to 
Snp = 0.1 the leading contribution is isovector in char¬ 
acter and the Lane parametrization works well. Even at 
S^p = 0.2 the second-order 5‘^p isoscalar contribution be¬ 
comes evident. At the highest value of 6np = 0.4, the 
isoscalar S‘^ term gives a contribution to the optical po- 













































































6 



FIG. 7: Isospin asymmetry dependence of the proton and 
neutron real optical potentials in infinite nuclear matter at 
saturation density po- The squares and circles are the cal¬ 
culated results, while the dashed and solid lines are lin¬ 
ear and quadratic hts, respectively. A scattering energy of 
E' = 50 MeV has been chosen. 


tential that is about 1/4 that of the isovector contribu¬ 
tion. 

The present work lays the foundation for improved 
modeling of nucleon-nucleus scattering away from the 
valley of stability. A strong isovector component to the 
imaginary part of the optical potential, as found in the 
present study, can inhibit radiative neutron-capture cross 
sections on exotic nuclei [60] and strongly influence r- 
process nucleosynthesis in cooler environments such as 
the tidally ejected matter in neutron star mergers. The 
isovector real part of the single-particle potential (to¬ 
gether with the isoscalar terms important at low 
energies) may be more relevant for neutron star inner 
crusts, where a lattice (or pasta structures) of neutron- 
rich nuclei interact with a background of free neutrons. 
At the higher energies attained at next-generation ra¬ 
dioactive beam facilities, our predicted sign change in 
the isovector real optical potential in the energy range 
110 < Einy < 200 MeV as well as the strength of the 
isovector imaginary optical potential Wj 8 — 12 MeV 
can be tested. 

Work supported in part by US DOE Grant No. DE- 
FG02-97ER-41014, the BMBE, the DEG cluster of excel¬ 
lence Origin and Structure of the Universe, by the DEG, 
NSEC (CRCllO). 


I. APPENDIX: 5lp CONTRIBUTIONS FROM 
THREE-NUCLEON FORCES 


short-distance contact term contribution, shown in Fig. 
3(a), has the form 


Uii{q,kf) = 


csk^f 


f 


127r4pAy 


( 20 ) 


The two diagrams from the one-pion exchange three- 
body force, Figs. 3(b) and 3(c), yield 


18(27r/^)4A,^ I 4u3 


■X 

3 u X x~^ 


u — X — X 


-1 


U l-^{u-\rxY 1 -\- {u — xY 

y (21) 


3 ^ 1 + (rt + x) 
+ —In ^ ^ 


2x 1 {u — xY 

The Hartree diagrams, Figs. 3(d) and 3(e), from the two- 
pion exchange three-body force proportional to ci and C 3 
give 




Su{2ci - C3) 

1 + 

1/1 A 2\ ^ / \ I ^ T (u, “h xY 

X ln(l + 4 uY + -(ci - C3) In -^ 

X 1 [u — xY 

Su{cs — Ci)(l + — u^) 

^ [1 + (rt + xY] [lE {u- xY] 

16u^{cs — 2ci)(l + 14 ^ — x‘^) I 

^[l + (u + a;) 2 ] 2 [l + (u-a;) 2]2 J' 


( 22 ) 


The Fock diagrams, Figs. 3(f) and 3(g), are split into two 
parts depending on ci and 03 ^ 4 : 


Uii{q,kf) = ^H{x, u) udlH{x,u) 

9 r u 

-2duH{x,u) + u[duH{x,u)] +dxH{u,x) -{3dx 
-2du)H{x,u)\^^^ -2H{u,u) + uH {u, u) dud^H {u, x) 

+ x) [u u) - 2duH{^, w)] I , (23) 


= ISiYfYx^ {|( 3 g 3 + 2 c 4 ) [duG{x, u)] ^ 

r ~\ 2u 

+{2c 4 - cs)G{x,u) 2duG{x,u) - udlG{x,u) + — (Scs 

J O 

-C4) [dul{x,u)]‘^ + 2(c3 + 04)1 {x,u) udll{x,u) 


Here we present results for the quadratic S^p correc¬ 
tions to the nucleon-nucleus optical potential from three- 
body forces at the Hartree-Fock level, which are found to 
be significantly smaller than the leading Snp terms. The 


1 / 7 / 

-2dul{x,u)^ +(c 3 - 2 c 4 ) d:cG{u,x)y-{3da;G{x,u)\^^ 
-2duG{x,u)\^^_^ - 2G{u,u)J + uG{u,u) dudxG{u,x) 
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+ J d^dxG(^,x)[udlG{i,u)-2duG{^,u)^ + 2(c3 

+C 4 ) dJ{u,x){^{Zdx - 2du)I{x,u)\^^^ - 21{u,u)^ 
+ul{u,u)dudxl{u,x) + J d^dxI{C,x)(udll{^,u) 
-2dul{^,u)) |, (24) 


with auxiliary functions defined in Eqs. (14)-(16). 
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